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ABSTRACT 


The interesting possibility of measuring the masses of high energy cosmic 
ray particles by observing pairs of extensive air showers arriving at the earth 
nearly simultaneously was proposed some years ago by Zatsepin (1951) and 
Gerasimova and Zatsepin (1960). Such showers would be created by the nuclear 
fragments originating as a result of the photodisintegration of massive nuclei 
interacting with the solar radiation held. In this paper we re-visit this possibility 
in the context of existing and proposed detectors of high and ultra-high energy 
cosmic rays considering a simple, yet realistic, model of the interplanetary 
magnetic held. The possibility of observing the mass fragmentation of cosmic 
rays directly, however, remains challenging. 

Subject headings: cosmic rays — interplanetary magnetic held — 
photodisintegration — extensive air showers 



1. Introduction 


Cosmic ray nuclei travelling towards the Earth can interact with the solar radiation 
field and photodisintegrate (Zatsepin, 1951), leading to the almost simultaneous creation 
of pairs of extensive air showers (EAS) in the earth’s atmosphere. If such pairs of showers 
could be detected and the initiating energies measured then the ratio of the greater energy 
to the lesser energy would give directly the mass of the heavier fragment. This assumes, as is 
most probable, that photodisintegration produces a single nucleon as one of the fragments. 
Thus in principle the mass of incoming cosmic rays could be measured rather directly. 
Gerasimova and Zatsepin (1960) studied this phenomenon in a simplified analytical way 
which was made necessary by the complexity of the problem and the limited computing 
power then available. As the problem of the mass composition of cosmic rays above the 
knee of the spectrum remains as important and controversial an issue as it was in the sixties 
(e.g., Watson 1997) it is of interest to assess the potential of any mechanism that could help 
to tackle the problem. With its remarkable simplicity, the Gerasimova-Zatsepin mechanism 
certainly falls in this category. 

Our main contribution is to address an important issue which was treated incompletely 
in the initial work namely the effect of the interplanetary magnetic field on the distribution 
of the expected core separations. We have made a significant improvement in this regard by 
performing exact orbit integrations in a realistic model of the interplanetary magnetic field. 
We also use modern estimates of the primary energy spectrum to calculate the expected 
rate of arrival of such pairs of showers and discuss the possibility of detection with current 
and planned instruments. 

The relevant photonuclear interactions are discussed, for example, in Danos and Fuller 
(1965), Hayward (1970) or Puget, Stecker and Bredekamp (1976). The energy range of 
interest in the rest frame of the cosmic ray nucleus spans from approximately 10 MeV up 
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to 150 MeV. However, it is in the region from ~ 10 MeV to ~ 30 MeV, the domain of the 
giant resonances of the nuclear photoeffect, where most of the interactions result in the 
emission of single nucleons, although two-nucleon emission can also occur but with much 
reduced probability. 

In the next section we describe our calculations and assumptions and discuss the 
results, while our brief concluding remarks are left to the last section. 


2. Numerical calculations and discnssion of results 

In dealing with this problem, two different aspects must be considered: the expected 
event rate and the relative deflection of the fragments. 

In the original paper (Gerasimova and Zatsepin 1960) a deviation of the order of 10“^ 
cm is quoted for an assumed homogeneous interplanetary field of the order of 10“® Gauss. 
This value is in error by several orders of magnitude for a typical particle with 7 ~ 10^ 
interacting with a photon in the vicinity of the Sun. An elementary calculation shows 
that deflections of hundreds of kilometers should be expected at Earth for a 7 ~ 10^ Fe 
nucleus losing a proton by photodisintegration at 1 AU. The oversight was subsequently 
recognized by Zatsepin, as reported by Ginzburg and Syrovatsky (1964), but no detailed 
calculations seem to have been made subsequently which include the effects of deflections 
in the interplanetary magnetic held. 

For our calculations we use a realistic, yet simple, model of the interplanetary magnetic 
held (IPMF) due to Akasofu, Gray and Lee (1980). This model takes into account four 
diherent components to describe the spiral 3-dimensional structure of the IPMF inside the 
central 20 AU of the Solar System: (1) the dipole component (B^jjpQ^g), (2) the sunspot 
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component (Bg^jjgpot); (3) the dynamo component and the ring current 

component (Bj.'j^g ). The total held is thus expressed as: 


^ “^dipole -^sunspot + dynamo -®ring (^) 

The solar dipole moment is well known (~ 3.4 x 10“^^ Gauss cm^). The dynamo 
component originates in a poloidal current system which exits the Sun at its poles, reaches 
the heliosphere at high latitudes, hows over this surface towards the ecliptic plane and 
hnally closes the circuit through an inward equatorial current. This system is supposed to 
be generated by a dynamo process induced by the solar rotation in the dipolar held. The 
ring current component is produced by a thin equatorial sheet current that extends up 
to the heliopause. The sunspot component constitutes the magnetic arcade immediately 
above the photosphere. In the model of Akasofu, Gray and Lee (1980) this component is 
represented by an ensemble of spherical dipoles just below the surface. Its main purpose 
is to allow the connection of all the held lines of the equatorial sheet to the Sun’s surface 
without signihcantly distorting the solar dipole held. Both the solar dipolar component 
and the sunspot component decrease as oc r~^ and therefore, outside the coronal region, 
the IPMF is mainly determined by the dynamo and ring current contributions. In hgure 
la we show the magnitudes of the total IPMF, and of the ring and dynamo components 
as a function of the cylindrical coordinate p over a plane located at z=0.01 AU above the 
ecliptic. Figure lb shows the corresponding cylindrical components of the IPMF. The 
signs over the curves indicate the sign of each component. Note, however, that all the 
components of the magnetic held reverse their direction with the 11 year solar cycle. 

The interplanetary medium is permeated by this magnetic held and by the photon 
radiation held originating at the Sun’s photosphere. For our purpose it is sufficient to 
assume that the radiation held spectrum is that of a black body at T^ff = 5770 K. 
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Consequently, the number density of photons at a distance tau (in astronomical units) from 
the Sun is: 


1 rtf 

~ ^ ^ ^ expiZ/al) - 1 

where e^v is the energy of the photons in eV in the reference system of the Sun. 

Every cosmic ray particle penetrating the Solar System must traverse this environment 
before reaching the Earth. In the rest frame of the nucleus the energy of the photons is 
boosted to: 




e 7 + V 7^ — 1 cos a ~ 2'ye cos — 


a 


( 3 ) 


where a is the angle between the momenta of photon and nucleus in the Sun’s reference 
system (e.g. Gerasimova and Zatsepin 1960). When the energy of the photons in the frame 
of the nuclei is larger than some few MeV, the nuclei can undergo photodisintegration. 
This process is most important between 15 and 25 MeV, in the region of the peak of the 
giant dipole resonance, although there is still a signihcant contribution to the cross section 
from energies beyond 25 MeV up to the threshold for photo-pion production at ~ 145 MeV 
(Puget, Stecker and Bredekamp 1976, Hillas, 1975). In the lower energy band, i.e. in the 
region of the giant dipole resonance (say, e<30 MeV), mainly single nucleons are emitted 
although double nucleon emission also takes place. At higher energies non-resonance 
processes are responsible for multinucleon emission. The corresponding cross-section 
parameters and branching ratios can be found in Puget, Stecker and Bredekamp (1976). 
Gerasimova and Zatsepin (1960) use a form for the cross-section that includes only 
the giant dipole resonance. We use the Gerasimova and Zatsepin representation of the 
cross-section at low energies, but have included additionally the non-resonant contribution 
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at higher energies up to the photo-pion production threshold. In Figure 2 we compare 
the cross-sections for Fe as calculated by Gerasimova an Zatsepin (1960), Puget, Stecker 
and Bredekamp (1976), an approximation by Hillas (1975), and the cross section used in 
the present work. It is evident that, although the Gerasimova and Zatsepin cross section 
is wider than the Lorentzian function used in Hillas’s approximation, it compares rather 
well with the combination of single and double nucleon emission as given by Puget and co¬ 
workers. Our approximation should give an upper limit to the photodisintegration rate. 
Uncertainties in the photodisintegration cross section will translate to only some few 10% 
and this does not alter our conclusion significantly. 

The emission of the nucleon(s) can be assumed to be isotropic in the reference system 
of the nucleus. Transforming to the Earth reference system, the emission of the fragments is 
concentrated within a cone of aperture ~ I /7 around the original direction of propagation 
of the parent nucleus. Therefore, at the high Lorentz factors of interest here (7 > 10^), we 
can assume that both fragments have, after the interaction with the photon, exactly the 
same direction as the incoming nucleus. Hence, we calculate relative deflections solely as 
the product of the differential bending of the fragments due to the action of the IPMF. This 
is, in fact, opposite to the approach taken by Gerasimova and Zatsepin. They neglected the 
effect of the IPMF and assumed that the distribution function of shower core separations 
was given only by the angular distribution of the fragments Lorentz-transformed to the 
Earth rest frame. 

We consider a spherical volume surrounding the Earth and extending up to r^ax = 4 
AU. The photon density is too low and the fragment deflections too large for any significant 
contribution to come from outside this region. A grid is constructed giving the fragment 
separation at Earth for a parent Fe nucleus interacting at any point inside the volume. As 
we are more interested in an upper limit than in an accurate calculation, we further assume 



that all of the cross-section goes into single nucleon emission producing, as daughters, both 
a Mn nucleus and a proton. 

To present the results we have chosen a polar coordinate system analogous to the 
galactic coordinate system. The reference system is centered on the Earth and its equatorial 
plane coincides with the ecliptic plane. The Sun is located at the origin of both longitude 0 
and latitude 9. Latitudes are positive to the North, while longitudes are positive to dusk 
and negative to dawn. Distances, in AU, are measured outwards from the Earth. 

In hgure 3 we show core separation distribution functions for three particular directions 
on the sky: noon (0 ~ 1.3 deg - i.e., a 5 R© perihelion), midnight (0 ~ —180 deg), 
and mid-afternoon (0 ~ 45 deg). For these plots the incident primaries are Fe nuclei 
of total energy E = 6.3 x 10^^ eV. This energy was chosen to be near the maximum of 
the fragmentation cross section for interaction with e ~ 1 eV photons (a typical energy 
of the solar radiation field photon). The three curves are proportional to the number of 
Gerasimova-Zatsepin events coming from each direction, and have been normalized such 
that the frequency is 1 for the smaller separation 6 arriving from the noon-side. It can 
be seen that the effect of the IPMF is much larger than that arising from the transverse 
separation of the fragments leaving the interaction; at 1 AU the angular spread of the 
fragments gives a separation of about 15 km. In fact the separation produced by the 
magnetic held is so large that the possibility of observing both partners of a disintegration 
process is rather limited at 7 ~ 10 ^ for any existing detector. 

The mean free path of a nucleus against photodisintegration is given by: 


1 

W) 



?r(/, e) X Cfrg 


2'ye cos^ 


«(/) 


X 2 cos^ 



(4) 


where I is the coordinate along the path of the nucleus, e is the photon energy in the rest 
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frame of the Sun, af^g is the fragmentation cross-section specihed previously, and a is the 
angle between the propagation directions of the nucleons and the photon, the latter being 
taken as a heliocentric radiovector. 

Therefore, dehning $00 as the unpertnrbed incoming cosmic ray flux at the external 
border, = 4A[/, and ^gz as the flnx of Gerasimova and Zatsepin fragment pairs, the 
relative GZ flnx is: 


where r^in ~ 
is carried ont. 


^G7 

Vgz = = 1 - exp 

rv~) 


dl 

W) 


( 5 ) 


0.02 AU is adopted as the inner spherical snrface np to which the integration 


Figures 4a and 4b show all sky maps of the ratio rjcz = ^gz/^od, i-e., of the fraction of 
Gerasimova-Zatsepin events among the incoming cosmic ray flux. The Sun is at the center 
of each hgure, and the ecliptic plane runs horizontally through the middle of the hgure 
(0 = 0). The shaded function in the backgronnd is f]GZ, while the contonr lines indicate 
the medians of the separation between the cores of correlated showers for each direction 
in the sky. Note that the labels on the separation contours of hgure 4a are logarithmic, 
while those in hgure 4b are linear. Both hgures are calculated for monoenergetic Fe nuclei 
at E = 6.3 X 10^^ eV (roughly the maximum of the photodisintegration cross section). In 
Figure 4a no acconnt has been taken of any acceptance ehects such as would be imposed in 
practice by a detector system: every event is connted irrespective of the separation of the 
showers. 


In Fignre 4b we show only those shower pairs for which the cores are separated by 
5 < 10 km. Snch a separation is relevant for the AGASA array (Ghiba et al. 1992) which 
has a collecting area of 100 km^. From these two hgures it can be clearly seen that high 
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values of rjcz only be obtained in the day-side in the vicinity of the Sun. Furthermore, large 
values of rjcz arise only when no cut-off in separation is considered: for any real instrument 
of hnite size (hgure 4b) the flux is several orders of magnitude smaller. In fact, so severe is 
this effect that when the separation is taken into account ( hgure 4b), the maximum rjcz is 
obtained on the night side of the earth around midnight. In other words, the small number 
of photons nightward of the earth is more than compensated by the high dehection on the 
day side which limits the effective integration volume to regions very near the Earth. The 
decrease of the effective integration volume on the day side can easily be seen by the change 
in symmetry of the function rjcz in going from hgure 4a to 4b. In the hrst case there is 
axial symmetry around the Sun because most of the events originates in its vicinity. In 
hgure 4b, when only small 6 events are accepted, they originate very near the Earth, and 
therefore there is symmetry with respect to the ecliptic plane, revealing the planar spiral 
topology of the IPMF in the neighborhood of the Earth. 

The low values of r] imply very low GZ huxes. To make an estimate we assume that 
all nuclei above E ~ 6 x 10^^ eV are of iron and that the integral hux is 3.8 x 10“^^ 
a value based on Fly’s Eye and Haverah Park data). Adopting an average 
value of riGz = 10“® (hgure 4b) we hnd the rate of GZ events to be ~ 0.01 per year on 
100 km'^ (the AGASA area). Even for the detectors of the Auger Observatory (3000 km'^ 
per site) only ~ 0.3 events per year with less than 10 km separation would land on the 
array. These rates are much too small for detection in the case of AGASA and for the Auger 
observatory the events will be difficult to identify within the background as the shower from 
the nucleonic fragment is likely to trigger only one detector because of the 1.5 km separation 
planned for the array. Our estimate of the rate is several orders of magnitude less than that 
made in the original examination of this ehect in part because of the magnetic dehections 
but also because modern estimates of the intensity at 6 x 10^^ eV are about 15 lower that 
believed in the late 1950s. 
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The situation does not get better at higher energies even for the planned large area 
detectors such as the Auger project and the proposed OWL satellites. The latter experiment 
involves a pair of twin satellites observing extensive showers in the upper atmosphere from 
the outer space. The advantage of such an experiment is the huge exposure area, 10^ — 10® 
km^. The disadvantage from the point of view of GZ events is, however, that OWL observes 
the night side of the atmosphere where the expected rjcz is rnuch smaller. Figures 5a and 5b 
illustrate the situation at higher energies: E = 3x 10^® eV and E = 3x 10^® eV respectively. 
In both cases a maximum separation 5max = 1000 km is used, representative of a 10® km^ 
experiment like OWL. The medians of the separations are conveniently smaller, but the GZ 
fluxes are again too low to be of practical use. Adopting an integral intensity of 2 x 10“^^ 
and rjcz = 10“^ gives only 0.06 events per year for an area on 10®A;m^. Hence 
the GZ flux is too low to make detection practical. Furthermore, at very high energies the 
separation between the cores is so small (5 ~ 5 km on the night side at 3 x 10^^ eV, or 
even 6 ~ 0.5 km ~ 10^® eV) that confusion might arise between the signals associated with 
showers. 


3. Conclusions 

We have re-analyzed the proposal of Gerasimova and Zatsepin (1960) of using pairs of 
correlated showers, originating in the photodisintegration of heavy nuclei interacting with 
the solar radiation held, as a mass-spectrometric technique. 

We consider a simple but realistic model of the interplanetary magnetic held to 
demonstrate that the magnetic dehection dominates the distribution function of core 
separations. The interactions with the solar radiation held inside a sphere of 4 AU around 
the Earth are calculated for incoming Fe nuclei under the simplifying assumption of single 



nucleon emission. The results are presented as all-sky maps, highlighting the considerable 
anisotropy of the solution. 

From our calculations it is apparent that the events arising from this very beautiful 
idea are too infrequent to be of use in any real experiment, either in operation or currently 
proposed, as a mass measuring technique. 


This work was done with the partial support of the Brazilian agency FAPESP. 



Figure Captions 


Figure 1: (a) Total magnitude of the interplanetary magnetic field and dynamo and 
ring current components as function of heliocentric distance over a plane parallel to the 
ecliptic plane and located at z=0.01 AU. (b) absolute value of the cylindrical components 
of the total magnetic held; the signs of the several components are indicated above the 
corresponding curves. 

Figure 2: Iron nuclei photodisintegration cross sections. See text for details. 

Figure 3: Core separation distribution functions. All the curves have the same 
normalization and can therefore be compared directly, each curve corresponds to a 
particular direction in the sky, as indicated in the respective labels, inside a solid angle of 
~ 2 X 10“^ sr. 

Figure 4: Fraction of GZ events, r]GZ, for Fe nuclei at E = 6.3 x 10^'^ eV (a) regardless 
of the separation 6 between correlated EAS and (b) only for 5 < 10 km. The contour lines 
indicate the median of the separation 6, and are labeled logarithmically (logio) in (a) and 
linearly in (b). The Sun is at the center of the image, while midnight is at the left and right 
borders of the hgures. 

Figure 5: As in hgure 4, but for (a) E = 3 x 10^® eV and (b) E = 3 x 10^® eV. The 
maximum separation allowed is in both cases Smax = 10^ km, and so this corresponds to a 
hypothetical experiment with an effective area of ~ 10® km^. The contour lines are labeled 
linearly with the median of 6 inside each solid angle. 
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